The Mott insulator — 10* order perturbation theory 
extended to infinite order using QMC 
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We present a new method, based on the combination of analytical and numerical techniques 
within the framework of the dynamical mean-field theory (DMFT) . Building upon numerically exact 
results obtained in an improved quantum Monte Carlo (QMC) scheme, lO"^ order strong-coupling 
perturbation theory for the Hubbard model on the Bethe lattice is extrapolated to infinite order. 
We obtain continuous estimates of energy E and double occupancy D with unprecedented precision 
0(10"^) for the Mott insulator above its stability edge Uci ~ 4.78 as well as critical exponents. The 
relevance for recent experiments on Cr-doped V2O3 is pointed out. 

PACS numbers: 71.10.Fd, 71.27.+a, 71.30.+h 



The exploration of the transition from a paramagnetic 
metallic phase to a paramagnetic Mott-Hubbard insula- 
tor is of fundamental interest in condensed matter theory 
Q . Essential features of this correlation phenomenon are 
captured by the single-band Hubbard model 



UD, 



(1) 



where the operator (cj^.) creates (annihilates) an elec- 
tron with spin a S {i, T} on site i and D = X^i'^iT^'U 
counts the number of doubly occupied sites; fiic — cj^c^^- 
Its nonperturbative solution in the thermodynamic limit 
is possible within the framework of the dynamical mean- 
field theory (DMFT) which maps the lattice problem 
onto an Anderson impurity model; this mapping becomes 
exact in the limit of infinite dimensionality Q. The na- 
ture of the Mott-Hubbard transition in the idealized case 
of a fully frustrated Hubbard model with semi-elliptic Q 
noninteracting density of states (DOS) has been a matter 
of intense debate As a result, a coexistence region 

of metallic and insulating solutions, enclosed by spinodal 
lines Uci{T) < Uc2{T), has been established below a crit- 
ical temperature T* « O.QGt. Only recently, estimates of 
energy E = {H)/N and double occupancy D = {D)/N 
per lattice site have been obtained in quantum Monte 
Carlo (QMC) calculations with high enough precision for 
reliably pinpointing the first-order metal-insulator transi- 
tion (MIT) line . A complete and universally accepted 
understanding of this prototype MIT is highly desir- 
able (i) as a prerequisite for controlled studies of more 
material-specific aspects of such transitions (e.g., band 
structure, degree of frustration) and (ii) since already the 
single-band description can explain surprisingly many as- 
pects of real-world experiments, which suggests a high 
degree of universality. In fact, recent experiments on 
V2O3 have not only found agreement with critical expo- 
nents at T* predicted by single-band DMFT but 



even observed well-defined spinodal lines. The interest- 
ing question to what extent these lines can be identified 
with those computed in purely electronic models is still 
open; one criterion will be determined in this Letter. 

In this work, we present results for the ground state 
insulating phase ("Mott insulator") with unprecedented 
precision. After introducing the Kato-Takahashi strong 
coupling perturbation theory (PT) used for an exact cal- 
culation of all contributions to the energy including order 
t^°/U^, we will sketch our advanced DMFT-QMC scheme 
which permits precisions of 0(10"^) (in E/t or D) at se- 
lected phase points. The error bounds will be confirmed 
at moderately strong coupling ([/ = 6t). Finally, the 
breakdown of finite-order PT observed at smaller U will 
be overcome using an extrapolation scheme to infinite 
order ( "ePT" ) which is controlled by its excellent agree- 
ment with QMC throughout the coexistence region. In 
addition to an estimate for Ud and to continuous esti- 
mates for E and D at all coupling strengths U > Ud, 
we will extract critical exponents from ePT (which are 
clearly out of reach for both QMC and plain PT). Beyond 
the direct physical interest, our results will provide an es- 
sential reference in the (now very active) development of 
new methods for solving (multi-band) DMFT, its cluster 
extensions, and hybrid schemes such as LDA-I-DMFT. 
The infinite-order extrapolation method ePT constructed 
in this Letter is similar in spirit to Singh's vastly suc- 
cessful extrapolation of lO"' order PT in J±/J\\ for the 
2d Heisenberg antiferromagnet to the isotropic limit ^3 ■ 
The common general ideas appear promising for many 
areas of physics, potentially including QCD. 

Perturbation theory - At half filling and in the absence 
of hopping, i.e. for tij = in QJ, all states without dou- 
ble occupancies have the minimal energy {U D) = 0; due 
to the arbitrariness of the spin direction on each site, the 
degeneracy of the space of unperturbed ground states 
is 2^. An expansion around this strong-coupling limit 
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FIG. 1: Low order diagrams in PT: A first hopping process 
creates an empty site (circle). After a total of 2n fioppings, 
impacting ?i other sites (spheres) , the ground state is restored. 



at T = u sing Kato's formalism was established b 
Takahashi for details in the case at hand, see Refs. 
and In this formalism, every ground state ipn of H is 
obtained by application of a well-defined operator T on 
a state (pn Uq: ipn = r<?5'n- Norm conservation implies 
T^T ~ Pq, where Pq is the projector on Uq. Thus, the 
Schrodinger equation HTjjn = E^ipn may be written as 



9n ■ 



h = ftf/f , 



where h acts as an effective operator in the unperturbed 
space and can be expanded in orders of t/U: 



^ — ^m— 1 



The operators hm generate electron transfers; since they 
have to preserve the absence of double occupancies, the 
associated diagrams are always closed. Specifically for 
the Bethe lattice, all closed paths are self-retracing. Nei- 
ther unlinked diagrams nor diagrams of odd order in t 
contribute (at half filling) for any lattice type 0. Fi- 
nally, graphs with more than two lines between any pair 
of sites are suppressed at least as XjZ. Fig. ^ illustrates 
the first three contributing orders. For the half-filled 



Hubbard model, 



order Kato- Takahashi PT yields 



coefficients for a generalized Heisenberg model includ- 
ing up to m-site interactions. In the paramagnetic phase 
and for Z = oo, an exact evaluation is possible since then 
all spin-spin correlation functions vanish. The lO"^ order 
result obtained from a computer implementation [isj 
of this diagrammatic approach will be discussed below. 

Quantum Monte Carlo - Deep enough within the in- 
sulating phase, all excitations of the insulator are expo- 
nentially suppressed. Consequently, ground state prop- 
erties can be estimated from measurements at finite (but 
small) temperatures T (and vice versa); the precision can 
be tested and further increased in explicit extrapolations 
T 0. At finite T, the impurity part of the DMFT self- 
consistency problem may be solved using the auxiliary 
field quantum Monte Carlo (QMC) method ^3; its ac- 
curacy is primarily limited by the discretization At of the 
imaginary time < t < T^^. In general, the At error 
consists of the Trotter error and an error associated with 
discrete Fourier transformations (FT) performed twice 
per iteration of the DMFT self-consistency cycle. The 
FT error is essentially eliminated in a new high-precision 
QMC algorithm where the low-frequency physics of the 
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FIG. 2; Ground state energy E and double occupancy D: 
PT according to and (jSJ (solid lines) in comparison with 
low-order PT (dashed/dotted lines) and QMC (circles). 



impurity is obtained from QMC while the high-frequency 
part is computed analytically. This method converges 
much faster for At — > at all frequencies than previ- 
ously used methods The QMC results shown below 
are based on simulations at T = 1/15 for discretizations 
0.09 < At < 0.25 (with quartic extrapolation At 0) 
using some 40 iterations with O(IO^) sweeps each. For 
U < 5.0, finitc-T corrections were estimated using addi- 
tional simulations at T = 1/20 and T = 1/25 [3. 

Results / - In the following, results are presented for 
the ground state (lower index omitted). PT for the 
half-filled Hubbard model Q with Bcthe DOS {W = At) 
yields the energy E{U) = Ept{U) + 0{t^^ C/""), where 
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for < = 1. Consequently, the double occupancy D{U) 
dE{U)/dU reads D{U) = Dpt{U) + 0{t^^ C/"^^), i.e.. 



Dpt{U) = 
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These PT results for E and D are shown as solid lines in 
the left and right panels of Fig. 13 respectively. The com- 
parisons with 2"^^ and 4**^ order PT (dashed/dotted lines) 
hint at a relatively fast convergence. In fact, a detailed 
analysis shows that the PT results at [/ = 6 (just beyond 
the coexistence region) are converged within about 10""" 
for both E and D. The agreement with QMC for this in- 
teraction confirms the correctness of the PT coefficients 
in ^ and of the QMC results within error bars. At 
smaller U, however, PT deviates markedly [by O(10~^) 
for D at U = 4.8] from the numerically exact QMC. In 
the remainder of this Letter, we will construct a new 
method, the extended perturbation theory (ePT), which 
yields continuous and thermodynamically consistent [l9l | 
estimates for E and D with even higher precision than 
QMC. In addition, we will obtain accurate estimates of 
the stability edge Ud and of critical exponents. 
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Extended perturbation theory ( ePT) - The idea of the 
ePT is an extrapolation of the coefficients a„ in the power 
series for the energy (for t = 1), 

1=1 

from a small number of known coefficients 02, ■ • • , a„Q to 
infinite order. Specializing; on the case where all a2i have 
the same sign |20|, we define the sequence 

Uci[2i + 2] = \/a2i+2/a2i- (5) 

If this sequence converges, the ratio criterion implies the 
convergence of series (0J for [/ > Ud = liuii^^ Ud 
and divergence below. More specificly, we assume 

Uci[n] = Uci + ui/n + U2/n^ + 0{n-^) . (6) 

Our assumption of this functional form will be strongly 
supported by fits to the known PT coefficients (cf. Fig.|2Jl. 
By virtue of such a fit, series (0} and its derivatives may 
be evaluated at arbitrary order. Comparisons of the nu- 
merical infinite-order limits of this scheme for E and D 
with QMC (cf. Fig. H and Fig. EJ will then give "nu- 
merical proofs" of the validity of Let us, however, 
first analyze the critical behavior: Eq. ^ and (jSJ imply 
Q!„ = a„C/J^" = n~'^[c + 0(1/??.)], where r = —ui/Uci- 
Re-expanding |@J around the critical point Ud , we obtain 

E{U) = Ud , «2. (1 + ; X ^ ^ , 

^1=1 Ud 

00 

= Ud [{ e„x") + e^_ia;^-^ -|- e^x^ + ...].(?) 

Tl = 

Note that any finite order in Q generates terms to all in- 
teger orders in x in (jT)). The leading singular contribution 
to E is associated with Cr-i = |r(l — r). For the numer- 
ical step within the ePT we use weighted [Ul quadratic 
least squares fits to © with either no adjustable param- 
eters ("unrestricted" fit) or with minimal tuning of Ud 
for enforcing analytic constraints or for obtaining error 
bounds from comparison with QMC. As seen in Fig. |S1 
the known coefficients (squares) are well represented by 
a "restricted" fit (solid fine) to be discussed below. Its 
small curvature as well as slight deviations from the PT 
values become apparent only in the inset, where an aver- 
age linear term has been subtracted. Also shown in the 
inset is the unrestricted fit (dashed line) which predicts 
Ud = 4.75 and r = 3.44. The shaded region marks ex- 
trapolation curves which lead to estimates for E and D 
compatible with QMC; thus, QMC implies error bounds 
4.72 < Ud < 4.79 and 3.36 < r < 3.53. Half-integer ex- 
ponents occur in the related Falicov-Kimball model 
and are generally expected in the context of mean-field 
theories. For our final estimates, we therefore assume 
an exponent r = 3.5 in a restricted fit which implies 




0.05 0.1 0.15 0.2 0.25 0.3 

1/n 

FIG. 3: Construction of ePT: PT values for Uc[n] (squares) 
are extrapolated to 1/n using a quadratic fit (solid line) 
restricted to t = 3.5. Values at smaller 1/n (circles) define 
ePT coefficients to all orders. Inset: restricted and unre- 
stricted fit (solid/dashed line), PT "data", and range com- 
patible with QMC (shaded) after subtraction of linear term. 
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FIG. 4: Main panel: Energy from QMC (circles), ePT (solid 
line), and truncated ePT (dashed/dotted lines); in all cases, 
the PT result [Eq. and solid line in Fig. |5] has been sub- 
tracted. Inset: energy correction in ePT beyond non-critical 
terms (solid line) in comparison with the leading critical term. 

Ud = 4.782 (in agreement with NRG HI and DMRG 
I2JI) and corresponds to the solid lines in Fig. 151 [25|. 

Results II - The ePT estimates for energy E and dou- 
ble occupancy D are shown in Fig. 0| and Fig. |S1 respec- 
tively (solid lines); in the main panels, the 10'^ order PT 
results have been subtracted for enhanced clarity. These 
results agree perfectly with QMC (circles) which con- 
firms the thermodynamic consistency of the latter [l^ . 
In contrast, low order truncated ePT (dotted/dashed 
lines) deviates markedly, in particular for D. The in- 
sets of Fig. 0] and Fig. [S] show contributions arising from 
the leading singular term in Q using 62.5 = —Cj^^/tt^ 
where c « —0.00725 is derived from the asymptotics of 
ri'^'^Q;[n]. These semi-analytic results agree well with the 
QMC/ePT estimates for contributions including and be- 
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FIG. 5: Double occupancy from QMC (circles) and ePT (solid 
lines) in complete analogy with Fig. 2] 

yond the leading critical term. Thus, the critical behavior 
is confirmed a posteriori by QMC. The polynomial (in x) 
representation of the ePT results used in the above 
analysis is more accurate than 10*'^ order PT for t/ < 6.6 
and should be used in future comparisons. 

Summary - We have combined the strengths of high- 
order FT and of an improved QMC scheme and have con- 
structed the ePT, an extension of PT to infinite order, in 
order to compute ground state properties of the Mott in- 
sulator within DMFT with unprecedented accuracy. We 
have confirmed earlier and independent estimates |23,|2J| 
for Uci, the lower (meta)stability edge of the insulator, 
and for the first time determined critical exponents (5/2 
for E, 3/2 for D). As both the spinodal line Uci{T) 
and the thermodynamic observables within the insulat- 
ing phase are nearly T-independent (for T < T*/2), these 
results are also excellent approximations at finite tem- 
peratures where the MIT and both spinodal lines have 
been observed experimentally in Cr doped V2O3 
In particular, a comparison between experimental esti- 
mates for the critical exponents and our predictions could 
shed light on the still ill-understood role of the lattice 
at the transition. With the authoritative results for E 
and D provided by ePT and QMC 113, aU details of the 
first-order Mott transition can now be clearly resolved 
which eliminates any room for alternative continuous- 
transition scenarios [3. Complementary results for the 
metallic ground state can be found in Ref. various 
finite-temperature results will be published elsewhere. 
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